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o 

We list scalar leptoquarks that mediate proton decay via renormalizable couplings to the Standard 



Model fermions. We employ a general basis of baryon number violating operators to parameterize 



contributions of each leptoquark towards proton decay. This then sets the stage for investigation of 
bounds on the leptoquark couplings to fermions with respect to the most current Super Kamiokande 
results on proton stability. We quantify if, and when, it is necessary to have leptoquark masses close 
to a scale of grand unification in the realistic SU(5) and flipped SU(b) frameworks. The most and 

Oh 

the least conservative lower bounds on the leptoquark masses are then presented. We furthermore 

single out a leptoquark without phenomenologically dangerous tree-level exchanges that might ex- 

plain discrepancy of the forward-backward asymmetries in tt production observed at Tevatron, if 

04 relatively light. The same state could also play significant role in explaining muon anomalous mag- 

> 

netic moment. We identify contributions of this leptoquark to dimension-six operators, mediated 
through a box diagram, and tree-level dimension-nine operators, that would destabilize proton if 
sizable leptoquark and diquark couplings were to be simultaneously present. 
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I. PROTON DECAY LEPTOQUARKS 

There has been a plethora of low-energy experiments capable of leptoquark discovery thus far. These 
have generated ever more stringent constraints on available parameter space for their existence. See, for 
example, [1-6] for some of the latest results. There also exists a large number of phenomenological studies 
of leptoquark signatures prompted primarily by various effects they could generate in flavor physics [7-1 1]. 
We are interested in a particular subset of scalar leptoquark states that are associated with proton decay. It 
is well-known that there exists only a small number of these states that can simultaneously violate baryon 
(B) and lepton (L) numbers [12-14]. The number of scalar leptoquarks that can mediate proton decay at 
the tree-level is even smaller [15]. Our aim is to present a comprehensive classification of leptoquarks and 
address a role these have in proton decay processes. 

Scalar leptoquarks that mediate proton decay certainly represent qualitatively new physics. Although the 
relevant operators associated with exchange of these states can be studied from an effective theory point of 
view, we prefer to trace their origins to a particular unification scenario in order to expose their dependence 
on underlying couplings. In fact we will study these states in two different unification frameworks that 
correspond to the SU(5) [16] and the flipped SU(5) [17-19], i.e., SU(5) x U(l), embeddings of the mater 
fields. These two scenarios are general enough to cover other possible embedding schemes. 

Let us start by spelling out qualitative differences between the scalar and vector, i.e., gauge boson, lepto- 
quarks that mediate proton decay at the so-called dimension-six (d = 6) level. (The latter have been studied 
much more extensively in the literature. See, for example, [12-14, 17, 20-24].) Firstly, vector leptoquarks 
comprise twenty-four states whereas the scalar ones comprise eighteen (fifteen) states in case neutrinos 
are Dirac (Majorana) particles. Secondly, whereas SU(5) contains only a half of all vector leptoquarks, 
the other half being in flipped SU (5), one can already find all possible proton-decay mediating scalar lep- 
toquarks in either SU (5) or flipped SU (5) framework. Hence, the scalar sector, although smaller, can 
potentially yield much richer structure with respect to the gauge one. Thirdly, the uncertainty in predictions 
for partial nucleon decay rates due to the gauge boson exchange resides entirely in a freedom to choose 
particular unitary rotations that need to be in agreement with observed mixing parameters in the fermionic 
sector as gauge bosons couple to matter with the gauge coupling strength. Scalar fields, on the other hand, 
couple to matter through Yukawa couplings. This brings additional uncertainties to potential predictions for 
relevant decay rates. 

The leptoquark states that simultaneously violate B and L quantum numbers tend to mediate proton de- 
cay at tree-level and are therefore taken to be very massive. However, we have investigated an SU(5) grand 
unified theory scenario [25] which resulted in a setup with a set of light leptoquarks. Namely, motivated 
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by the need to explain anomalous events in ti production at Tevatron [26, 27], we have found that a light 
color triplet weak singlet scalar could contribute to ti production and explain the observed increase of the 
forward-backward asymmetry [28]. We have accordingly demonstrated that the unification of the funda- 
mental interactions is possible if that set of light scalars is a part of the 45-dimensional representation [28]. 

In flavor physics, due to recent accurate measurements at Tevatron and LHCb, the presence of new 
physics (NP) in B systems seems rather unlikely. The muon anomalous magnetic moment, on the other 
hand, still leaves some room for NP contributions. The impact of potentially light leptoquark scalars, 
including the light color triplet weak singlet scalar, on the low energy and hadron collider phenomenology 
within that context has been investigated in Refs. [28-32]. 

The color triplet weak singlet scalar state we have singled out does not generate proton decay at the tree- 
level. However, one can still construct, as we show later, higher order loop diagrams that yield effective 
d = 6 and tree-level d = 9 operators which can destabilize proton. The natural question then is whether 
one can simultaneously address the ti asymmetry and the muon anomalous magnetic moment by using the 
very same leptoquark. We investigate this issue in detail in Sec. VI. 

This paper is organized as follows. In Sections II and III we list all proton decay inducing leptoquarks in 
SU(5) and flipped 577(5) unification frameworks and specify their Yukawa couplings to the SM fermions. 
In Sec. IV we introduce the effective dimension-six operators for proton decay and calculate associated 
effective coefficients for each leptoquark state. Sec. V is devoted to a study of conservative lower bounds on 
the color triplet leptoquark mass within phenomenologically realistic 577(5) and flipped SU(5) scenarios. 
In Sec. VI we study leptoquarks that do not contribute to proton decay operators of dimension-6 at tree-level. 
We conclude in Sec. VII. 

II. LEPTOQUARKS IN SU(5) 

The scalars that couple to matter at tree-level reside in the 5-, 10-, 15-, 45- and 50-dimensional rep- 
resentations of SU(5) because the SM matter fields comprise 10; and 5j, where i,j = 1,2,3 repre- 
sent family indices. Namely, 10; = (1, 1, 1), © (3, 1, -2/3); © (3, 2, 1/6); = (ef,uf, Q { ) and 5j = 
(1, 2, -l/2)j © (3, 1, 1/3) j = (Lj,df), where Qi = (u; di) T and Lj = (vj ej) T [16]. Possible con- 
tractions of the matter field representations hence read 10 (g> 10 = 5 © 45 © 50, 10 © 5 = 5 © 45 and 
5 © 5 = 10 © 15. Theory also allows for addition of right-handed neutrinos that can be introduced, for 
example, in the form of SU(5) fermionic singlets (1) without the need to enlarge the scalar sector. Note that 
one can also introduce additional non-trivial representations of matter to generate observed fermion mass 
parameters in the lepton [33] and quark [34] sectors. That, however, would not alter our operator analysis 
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for large enough masses of extra matter fields. 

Relevant decomposition of scalar representations to the SM gauge group, i.e., 517(5) — > SU (3) x 
SU(2) x 17(1), is given below [35]: 

• 5 = (1,2,1/2)0(3,1,-1/3); 

• 10 = (1, 1, 1) (3, 1, -2/3) (3, 2, 1/6); 

• 15 = (1, 3, 1) © (3, 2, 1/6) © (6, 1, -2/3); 

• 45 = (8, 2, 1/2) © (6, 1, -1/3) © (3, 3, -1/3) © (3, 2, -7/6) © (3, 1, -1/3) © (3, 1, 4/3) © 
(1,2,1/2); 

• 50 = (8, 2, 1/2) © (6, 1, 4/3) © (6, 3, -1/3) © (3, 2, -7/6) © (3, 1, -1/3) © (1, 1, -2). 

Only 5, 15 and 45 contain electrically neutral components and are thus capable of developing phe- 
nomenologically viable vacuum expectation values (VEVs). Contributions to the up-quark, down-quark 
and charged lepton masses can come from both 5 and 45 whereas Majorana (Dirac) masses for neutrinos 
can be generated by VEV of 15 (5). 

The scalar leptoquark states that violate both B and L quantum numbers are (3, 1, —1/3), (3, 3, —1/3) 
and (3, 1,4/3), if one assumes neutrinos to be Majorana particles. These states reside in 5, 45 and 50. 
However, if one allows for the possibility that neutrinos are Dirac particles there is another leptoquark — 
(3, 1, —2/3) — that is found in the 10 of SU(5) that violates both B and L and could thus also destabilize 
proton. To that end we consider both the Majorana and Dirac neutrino cases to keep the analysis as general 
as possible. Altogether, there are eighteen (fifteen) scalar leptoquarks that could mediate proton decay 
in case neutrinos are Dirac (Majorana) particles. The leptoquarks in question are all triplets of color as 
they must contract with lepton and quark states into an SU(3) singlet. Tables I, II, III and IV summarize 
couplings to the matter of relevant states that reside in 50-, 45-, 10- and 5-dimensional representations, 
respectively. 

We observe that in the SU(5) framework the primary obstacle to the proton stability seems to be the need 
to generate Yukawa couplings relevant for the charged lepton and down-quark masses. These receive equally 
important contributions from the 10/5^ 5* and 10^ 45* contractions [36]. It is clear from Tables II and IV 
that both of these, individually, generate potentially dangerous couplings. The up-quark Yukawa coupling 
generation, on the other hand, via the lOjlO.,45 operator seems not to pose any danger whatsoever as can 
be seen from the second column in Table II. However, that operator cannot generate viable masses for all up- 
quarks due to the antisymmetry of the corresponding mass matrix. The lOjlO-,5 contraction does provide 
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SU(5) 


YtflOi 10,50 


(3,1,-1/3) 
A 


12-V2 eobcK io + ^io ]d T. Cub . Ac 

3-1/2^10 + ^10 ]e CT C7u C Aa 



TABLE I. Yukawa couplings of the B and L violating scalar in the 50-dimensional representation of SU(5). 
a,b,c — 1,2,3 = 1,2,3) are color (flavor) indices. are Yukawa matrix elements associated with the 

relevant contraction in the group space of SU(b). 



SU(8) 


y^iOiio^s 


10,5,45* 


(3,1,-1/3) 
A 


2^ K io_^o ]ef T Cn c Ao 


2 1 Y^e abc u a - r Cd b: jA* 
-2-^f^C ei A* 


(3,3,-1/3) 
(A 1 , A 2 , A 3 ) 


2^e abc [Y™-Yj?]dZ i Cd bj Al 




-2e abc [Y±° - Y £]d T ai Cu bj Al 


2-V2^T CejAr 


-2i/\ abc [Ytf -YV>}u T ai Cu bj A* c 


-yfCCe.Af 


(3,1,4/3) 
A 







TABLE II. Yukawa couplings of the B and L violating scalars in the 45-dimensional representation of SU(5). a,b,c = 
1, 2, 3 (i, j — 1, 2, 3) are color (flavor) indices. Y^° and Y^ are Yukawa matrix elements. 

viable up-quark masses but the price to pay is resurrection of the proton decay issue. To conclude, the only 
operator that can be considered innocuous in the Majorana neutrino case is the 10j 10^45 contraction. 

III. LEPTOQUARKS IN FLIPPED SU(5) 

Another possibility to unify the SM matter into an 5t7(5)-based framework leads to the so-called flipped 
SU(5) scenario [17-19]. A single family of matter fields in flipped 577(5) can be seen as originating from a 
16-dimensional representation of 50(10). Actually, flipped SU(5) is not necessarily completely embedded 
into 5O(10). Nevertheless, the generator of electric charge in flipped SU(5) is given as a linear combination 
of a f7(l) generator that resides in SU(5) and an extra J7(l) generator as if both of these originate from an 
50(10) — > 577(5) x U(l) decomposition. This guarantees anomaly cancelation at the price of introducing 
one extra state per family, i.e., the right-handed neutrino v c . The transition between the 577(5) and flipped 



6 



SU(5) 


yAiOiij-io* 


^5,5,10 


(3,1,-2/3) 
A 




2- 1 / 2 ea6c ^f^C^.A c 

y5 _ „y5T 



TABLE III. Yukawa couplings of the B and L violating scalar in the 10-dimensional representation of SU(5). a,b,c = 
1, 2, 3 (i, j = 1, 2, 3) are color (flavor) indices. Yh and 5^!j are Yukawa matrix elements. 



SU(5) 




if 10^- 5* 


^■5,1,5 


(3,1,-1/3) 
A 


2e abc [Y z f +Y^]d^Cu bj A c 
-2[Y$ + Y$]e<? T Cu c aj A a 


-2-V^|^^,A* 





TABLE IV. Yukawa couplings of the £> and L violating scalar in the 5-dimensional representation of SU (5). a,b,c = 
1,2,3 (i, j = 1,2, 3) are color (flavor) indices. Y^ , and Yl are Yukawa matrix elements. 

SU(5) embeddings is then provided by d c i/ 7 , e c i/ , u ^ d and i> •<-)■ e transformations. Flipped 
SU(5) thus predicts existence of three right-handed neutrinos as these transform nontrivially under the 
underlying gauge symmetry. 

The matter fields in flipped SU (5) comprise 10+ l , 5 { 3 and where the superscripts correspond to 
the extra U(l) charge assignment. To obtain the SM hypercharge Y one uses the relation Y = (Y(U(1)) — 
Y(U(l)su(5)))/5, where Y(U(1)) and Y(U(l)su(5)) represent the quantum numbers of the extra U(l) 
and the U (1) in SU(6)(-¥ 577(3) x SU (2) x U(l)), respectively. One then obtains the electric charge as 
usual, Q = Y + T 3 . 

The scalar sector that can couple to matter directly is made out of 50 2 , 45 -2 , 15 +6 , 10 +6 , 5~ 2 and 
1~ 10 . Representations that can generate contributions to the charged fermion masses and Dirac neutrino 
masses are 45 2 and 5~ 2 , whereas Majorana mass for neutrinos can originate from interactions with 15 +6 . 
Leptoquarks that violate B and L reside in 50 2 , 45~ 2 , 10 +6 and 5~ 2 with relevant couplings to matter 
given in Tables V, VI, VIII and VII, respectively. 

In flipped SU(5) the main obstacle to matter stability is the the generation of the up-quark masses. 
Namely, these can be generated through 10+ ^ 3 5* +2 and/or 10f 1 5j 3 45* +2 contractions. Both of these 
are dangerous as far as the proton decay is concerned as can be seen from Tables VI and VIII. All other 
contractions, in the Majorana neutrino case, are actually innocuous. 
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SU(5) x 17(1) 


YAOlO+ilo+^O" 2 


(3,l,-l/3)- 2 
A 


12-V\ abc [YV + YV]u T ai Cd bj A c 
S-^[Y^ + Y^rCd^A a 



TABLE V. Yukawa couplings of the B and L violating scalar in 50-dimensional representation of SU(5). a, b, c = 
1, 2, 3 (i, j — 1, 2, 3) are color (flavor) indices. y7° are Yukawa matrix elements. 



SU(5) x (7(1) 


y^io+'io+Ub- 2 


y/}l(L5/ 3 45*+ 2 


(3,l,-l/3)- 2 
A 


2^[Y l f-Y^]vf T Cd c aj A a 


2-^e abc d c JCu^A* c 
2- 1 Yf J u T ai Ce J Al 


(3,3,-l/3)- 2 
(A 1 , A 2 , A 3 ) 


2 1 ' 2 e abc [Y l f - Y^]u T al Cu bj Al 


Y^CejAl* 


-2e abc [Y±° - Y™]u T al Cd b3 Al 


2-^Y^CvjAl* 
2-^Y^ at Ce,Al* 


-2^e abc [Y^-Y^]d T al Cd bj Al 




(3,l,4/3)- 2 
A 


2^[Y l f^Y^]e abc df a T Cd^A c 


-Y^rCu^K 



TABLE VI. Yukawa couplings of the B and L violating scalars in 45-dimensional representation of SU(5). a,b,c = 
1, 2, 3 (i, j — 1, 2, 3) are color (flavor) indices. Y±° and are Yukawa matrix elements. 



IV. PROTON DECAY 



Let us discuss proton decay operators due to the scalar leptoquark exchange of the lowest possible 
dimension in detail. These are dimension-six operators made out of three quarks and a lepton that violate 
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SU(5) x (7(1) 


FAlO+ 1 l+ 5 10*- 6 


Y$S t % ' lO " 


(3,l,-2/3)+ 6 
A 




y5 — _y5T 



TABLE VII. Yukawa couplings of the B and L violating scalar in 10-dimensional representation of SU(5). a, b, c 
1, 2, 3 (i, j — 1, 2, 3) are color (flavor) indices. yl and y^ are Yukawa matrix elements. 



SU(5) x £7(1) 


FA°10+ 1 10+ 1 5- 2 


y i fio+ 1 5" 3 5*+ 2 




(3,l,-l/3)- 2 
A 


-2e ahc [Y™ +Y™]ul i Cd hj A c 


-2-^Y^Ce.K 


yAu^CefA a 



TABLE VIII. Yukawa couplings of the £? and L violating scalar in 5-dimensional representation of SU(5). a, b, c 
1, 2, 3 (z, i = 1, 2, 3) are color (flavor) indices. Y±°, Yfi and y>- are Yukawa matrix elements. 



B and L by 1 unit. They are summarized below 



OH(d a ,e/s) 
H (d a ,e%) 
H {d%,ep) 

H (d a ,di3,Vi) 
H (d a ,d^,Vi) 

H {d a ,d c ^v?) 

H {d c al d c ,,uf) 



a{d a ,ep) u T L C 1 d a u T L C 1 ep, 
a(d a , e c p ) u T L CT 1 d a ef L C^u * , 
a{d%, ep) d£ L C' 1 u c * u T L C~ x ep , 
a(d%, eg) df L C' 1 u c * ef L C^u * , 
a(d a , dp, Vi) u T L C^ 1 d a da L C~ x V{ , 
a(d a , dp, dp L C~ x u c * d T a L CT 1 v x , 
a(d a , dFp, v?) u T L C~ l d a vf L C' 1 d c p * 
a(d£, d%, v?) df L C~ x u c * vf L C 



- 1 d c * 



(1) 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



Here, i(= 1, 2, 3) and a, /3(= 1, 2) are generation indices, where all operators that involve a neutrino are 
bound to have a + (3 < 4 due to kinematical constraints. L{= (1 — 7s)/2) is the left projection operator. 
The SU (3) color indices are not shown since the antisymmetric contraction e a b c q a qbq c is common to all the 
above listed operators. This notation has already been introduced in Ref. [24]. 

These operators allow one to write down explicitly d = 6 proton decay contributions due to a particular 
leptoquark exchange [24]. To that end we first specify our convention for the redefinition of the fermion 
fields that yields the up-, down-quark and charged lepton mass matrices in physical basis: Mu,d,e 



M^ E . These are U£,M V U = M^, D L C M D D = M%*, and E£M E E = M^. The quark mixing is 



rdiag 



diag 



rdiag 
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U^D = Vud = K\VckmK2, where K\ and K2 are diagonal matrices containing three and two phases, 
respectively. In the neutrino sector we have N^M^N = M^ ag {N T M^N = M^ ag ) in the case of Dirac 
(Majorana) neutrinos. The leptonic mixing E'N = Ven = K^VpunsK^ in case of Dirac neutrino, or 
Ven = KsVpmnS i n the Majorana case. K3 is a diagonal matrix containing three phases whereas 
contains two phases. Vqk m (Vpmns) is the Cabibbo-Kobayashi-Maskawa (Pontecorvo-Maki-Nakagawa- 
Sakata) mixing matrix. 



A. Tree-level exchange (d = 6) operators in SU(5) 



The only relevant coefficient for A = (3, 1, —1/3) from 50 is 



a(d a , e 



(U T (Y 10 + Y w 1 )D) la (EUY W + Y 



6m 2 A 



(9) 



where toa is a mass of leptoquark in question. (See Table I for details on notation for Yukawa couplings of 
the 50-dimensional representation to the matter.) 

The relevant coefficients for A = (3, 1, —1/3) from 45 are 

1 



a{d^,e B ) 
a(d^,ef) 
a(d a ,dp,Vi) 



1 



771,/ 



{tf c Y^U* c U (U 1 Y E)i/3 , 
(rf c Y^U* c ) al (EUY l0 -Y 10T )^ c ) pi , 



{tf c Y^U* c )pi (D T Y b N) 



The relevant coefficients for A = (3, 1, —1/3) from 5 are 



a(d a ,e B 



a(d a , e B 



a(d^,e g 



a{d c a ,e c : 



a(d a ,d/3,Vi 
a(d Q ,d^,Ui 
a(d a ,d c B ,Vi 
a{d c a A C B ^ 



V ' 2 (U T (Y 10 + Y WT )D) la (U T Y"E) 1B 



irr 



m, 



(U T (Y W + Y^ )D) la {El{Y w + Y^yU* c ) B1 



-10 T 



t f V 10 



r 10 T\ f 7 



(rf c Y^U* c ) al (U r Y»E) 13 , 



2m\ 



in 



V2 



(U T (Y 10 + Y W1 )D) la {D 1 Y b N) B 



Tv5i 



77?- 



(rf c Y^U* c ) B1 (D r Y"N) ai , 



2m\ 



777, 



-(U T (Y w + Y WT )D) la {NlY^D* c ) iB , 



t yitn* 



(10) 

(11) 
(12) 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
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The relevant coefficients for A 2 G (3, 3,-1/3) from 45 are 



a(d a , ep) = (U T (Y W - Y WT )D) la {U T Y*E) W , (21 , 



a(d a , dp, Vi ) = ~iS-(U T (Y 10 - Y 10T )D) la (D T Y 5 N)/3i . (22) 
A 2 



The relevant coefficient for A 1 G (3, 3,-1/3) from 45 is 

a(da,d fi ,vi) = ^{U T Y^N) Xi (D T (Y 10 -Y WT )D)p a , (23) 

M A 1 

where the extra factor of 2 comes from two terms in Fierz transformation 

(s^Ld) {yCLu) = - (vPLs) (y°Ld) - (up ' Ld) (y^Ls) . (24) 
The only relevant coefficient for A = (3, 1, —2/3) from 10 is 



a{d c a1 d%vf) = -_ L^. {Dl(Y*-Y^D* c ) Pa {N^U^a . (25) 



Finally, A 3 G (3, 3, -1/3) and (3, 1, 4/3), both from 45 of SU(5), do not contribute to proton de- 
cay at tree-level. This is due to antisymmetry, in flavor space, of their couplings to the pair of up-quarks. 
Nevertheless, both states still induce proton decay through loops at an effective d = 6 level. We present 
systematic study of these contributions for the (3, 1,4/3) case in Section VI. There we also spell out con- 
tributions of the (3,1,4/3) leptoquark to dimension-nine tree-level proton decay amplitudes. Equivalent 
contributions of A 3 G (3, 3, —1/3) are not pursued since the components A 1 and A 2 from the same state 
already contribute at leading order. In this manner, higher order contributions of A 1 ' 2,3 would only play a 
role of radiative corrections. 

B. Tree-level exchange (d = 6) operators in flipped ST/ (5) 

The only relevant coefficient for A = (3, 1, — 1/3)~ 2 from 50 2 is 

a{d a ,d%vf) = JL (U T (Y w +Y WT )D) la (NUY 10 +Y WT ?D* o ) ip . (26) 
The relevant coefficients for A = (3, 1, — l/3)~ 2 from 45 -2 are 

a ( d a , e/j) = ^ (D^*U* c ) al {U T Y^E) W , (27) 

a(d a ,d%,vi) = (D j c Y~ 5 *U* c ) pi (D T Y^N) a i , (28) 

a(d£,$,«F) = -j2-^(D* c Y~ 5 *U* c )^ (iVj (F 10 - Y WT ^D* c ) ia . (29) 
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The relevant coefficients for A = (3, 1, —1/3) 2 from 5 2 are 

V2 



a{d a ,ep 

4 



a(d a , e c 



a{d a , ep 
a{d c a ,e c 



a(d a ,dp,Vi 
a(d a ,d C p,Vi 
a(d a ,dC,vf 

a(dC,d C p,v? 



(U T (Y w + Y 10T )D) la (U^E)^, 



TV5-, 



m. 



(U T (Y 10 + Y 10T )D) la (ElY^U* c )p 



t v 1 Ut*' 



2m 



, (D ] c Y^*U* c ) a i (U T Y 5 E) W 



V2m 2 A 



m 



^(U T (Y W + Y 10T )D) la (N T Y^ T D)ip . 
A 
1 



7V5; 



2771^ 
-4 



(t^yio + y 10T )z?)i Q (A£(y iU + y wt )W* c ) 



m, 



V2 



2 (z^y 5 *[/£) 1/3 (4(y 10 + y 10T )tz£) fa 



The relevant coefficients for A 2 G (3, 3, —1/3) 2 from 45 2 are 

V^2 



a(d a ,ep) 
a(d a ,dp,Ui) 



(U T (Y w -Y WT )D) la (U r Y*E) lf3 



T v 5; 



m a2 



7V5; 



The relevant coefficient for A 1 £ (3, 3, -l/3)~ 2 from 45~ 2 is 

2V2 



a(d a ,dp,Vi 



Mil 



3 (c/^^H^cy 10 -^^)^)^ 



where the extra factor of 2 comes from Fierz transformation (24). 
The only relevant coefficient for A = (3, 1, 4/3) ~ 2 from 45 -2 is 



a{dC,d c p ,vf 



(30) 
(3D 
(32) 
(33) 

(34) 
(35) 
(36) 

(37) 

(38) 
(39) 

(40) 



(41) 



where the extra factor of 2 again comes from Fierz transformation. 

The relevant coefficients for (3, 1, -2/3)+ 6 G 10 +6 and A 3 G (3, 3, -1/3)" 2 G 45~ 2 are not present 
at the leading order due to antisymmetry of the couplings to the up-quark pair. The higher order contribu- 
tions of the former state are discussed in Section VI. 



V. LEADING ORDER CONTRIBUTIONS 



Scalar fields couple to matter through Yukawa couplings. This introduces uncertainties to predictions 
related to any process that involves scalar leptoquark exchange. It is thus natural to ask if, and when, it is 
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necessary to have leptoquark masses close to a scale of grand unification. We will address this issue in the 
££7(5) and flipped SU (5) frameworks in what follows. 

A. Color triplets of 577(5) 

If the Yukawa sector relevant for proton decay through scalar exchange is not related to the origin of 
fermion masses and mixing parameters one cannot make any firm predictions. For example, all operators 
that correspond to the exchange of the triplet scalar in the 5-dimensional representation of ££7(5) can be 
completely suppressed if (U T (Y W + Y WT )D)i a = and {D ] C Y^^U* C ) al = 0, a = 1, 2, in the Majorana 
neutrino case. (Recall, it was the exchange of this scalar that has led to the so-called doublet-triplet splitting 
problem within the context of the Georgi-Glashow ££7(5) model [16].) The suppression is certainly viable 
if the entries of £7, Uc, D, Dq, Y 10 and Y 5 are all free parameters. The first set of conditions can be insured 
if, for example, Y 10 = — Y 10T . This solution has already been pointed out in Ref. [24]. The second set of 
conditions can also be easily satisfied although what we find defers from what has been presented in [24]. 

Things, however, change in models where the connection between Yukawa sector and fermion masses 
is strong. Let us thus analyze predictions of the simplest of all possible renormalizable models based on 
the SU (5) gauge symmetry. We want to find what the current experimental bounds on the partial proton 
lifetimes for processes presented in Table IX imply for the masses of color triplets if the theory is to be 
viable with regard to the fermion mass generation. We analyze all these decay modes to make our study as 
complete as possible. 



PROCESS 


t p (10 33 years) 


p — > 7r°e + 


13.0 [37] 


P — > 7T°/Z + 


11.0 [38] 


p -> K°e+ 


1.0 [39] 


p -> K°fj,+ 


1.3 [39] 


p — > r/e + 


4.2 [40] 


p — > rjfi + 


1.3 [40] 


p — > ir + i> 


0.025 [41] 


p — S- K + v 


4.0 [37] 



TABLE IX. Experimental bounds on selected partial proton decay lifetimes @ 90 % CL. 

We demand in what follows that the theory is renormalizable and thus neglect possibility that higher- 
dimensional terms contribute to (super)potential at any level. We furthermore take the simplest possibility 
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for the generation of phenomenologically viable fermion masses and mixing parameters. Namely, we de- 
mand that both 5 and 45 of Higgs contribute to the down-quark and charged lepton masses [36]. We further 
take all mass matrices to be symmetric, i.e., Mu,D,E = Mjj D E . This then allows us to consider two par- 
ticular scenarios. The first (second) one represents the case when the contributions of the (3, 1, —1/3) state 
from the 5-dimensional (45-dimensional) representation dominates. Our analysis is self-consistent as the 
symmetric mass matrix assumption eliminates contributions to proton decay of all other color triplets. Note 
also that any mixing between the triplets can be accounted for by simple rescaling of relevant operators. 



1. The charged anti-lepton final state 

We start our analysis with proton decay due to exchange of the triplet state from the 5-dimensional 
representation. To find widths for the charged anti-leptons in the final state one needs to determine a(d a , ep), 
a(d a , e^), a(da, ep) and a(d^, e^). If the Yukawa couplings are symmetric the relevant input for these 
coefficients reads 

{u t {y io + Y WT )D) la = (M$^V UD ) la , (42) 

V2v 5 

{U T Y*E) W = m D M A ^vl D U* 2 + U 2 M^) W , (43) 

(D^U*)ai = ~ (3Mj iag Vf D + Vl D USM%*ulU , (44) 

( £ t( ylO + yiOTjf^ = _1 {U T M ^ hl t (45) 

V2v 5 

where JJ% = U T E* and v$ represents the VEV of the 5-dimensional representation. U2 entries and 
V5 are primary sources of uncertainty. Our normalization is such that | ^5 1 2 / 2 + 12|v45| 2 = v 2 , where 
v(= 246 GeV) stands for the electroweak VEV. W45 is the VEV in the 45-dimensional representation. A 
connection between Yukawa couplings and charged fermion mass matrices is spelled out elsewhere [31]. 

We outline details of our calculation using the p — > e^ir channels. Here, S = 1 (5 = 2) corresponds 
to e + (// + ) in the final state. The decay widths formulas we use are summarized in Ref. [24]. For these 
particular channels we have 

r(p ejn ) = (Mdi, eg) + /3o(df , e 5 )\ 2 + \aa(d lt e c s ) + /3a(df , ef)\ 2 ) (1+D+F) 2 , 

(46) 

where a and /3 are the so-called nucleon matrix elements. (See Section VI for more details on a and /3.) 
F+D and F — D combinations are extracted from the nucleon axial charge and form factors in semileptonic 
hyperon decays, respectively [42, 43]. We take f n = 130 MeV, m p = 938.3MeV, D = 0.80(1), F = 
0.47(1) and a = -(3 = -0.0112(25) GeV 3 [43]. 
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The uncertainty in predicting partial decay rates persists even in the minimal SU (5) scenario with sym- 



metric Yukawa couplings. This is evident from the U 2 dependence of p — > e^ir partial decay widths 



T(p e+vr ) 



Op - "V 



2 \2 







(V UD ) u [m u + - A m d ] + |(^ D C/|M^""T^ ) , , 



T diag TT \\ 



\m D M^ D W 2 ) ls + \{U 2 M d ™ 9 ) u 



+ Mm u (U 2 ) 15 \ 2 (l + D + Ff. 



One can suppress (enhance) F(p — > e~$ir°) with regard to U 2 numerically to generate the least (most) con- 
servative lower bound on the mass of the scalar triplet in the 5-dimensional representation of SU (5). This, 
however, should be done simultaneously with all other partial decay modes to generate a self-consistent 
solution. We will do that after we address proton decay into neutral anti-leptons in the final state. The 
important point is that even in the case of symmetric Yukawa couplings one cannot test the SU(5) theory 
when the scalar triplet exchange dominates. This is in stark contrast to what one obtains for the gauge d = 6 
contributions in the SU(5) framework [20]. 

Nevertheless, we can already outline how one can find a maximum of F(p — > e^7r°) with regard to U 2 
to obtain the most conservative bound, from the model building point of view, on ttia without resorting to 
elaborate numerical analysis. The idea is to have Yukawa couplings of the third generation contribute as 
much as possible towards relevant amplitudes. The only possibility to achieve that is to make the 1 1 element 
of U 2 Mg ag U 2 matrix in Eq. (44) as large as possible. This can be done with a simple ansatz 

/ l\ 
10 



Uo = U 1 E* 



1 



(47) 



to obtain the following simplified expressions 

r(p 



e 7r 



r(p /i + vr° 



3 2 m p a 2 , 

2 l MMmi \(VUD. 

3 2 m p a 2 , 

2 l MKmi \(VUD. 
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12 



m 2 b m 2 T {l + D + Ff 



m 2 s m 2 T {l + D + Ff 



(48) 
(49) 



Interestingly enough, these expressions when combined with experimental input yield comparable bounds 
on the mass of the leptoquark in question. 

In the previous analysis it was assumed that contributions to proton decay of the triplet in the 5- 
dimensional representation dominate over contributions of triplets in the 45-dimensional representation. 
Let us now see if and when that is truly the case. Our assumption that the mass matrices are symmetric 
implies that the only other contribution to p — > e^7r° (5 = 1,2) channels originates from exchange of the 
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(3, 1, —1/3) state in the 45-dimensional representation. In fact, the only relevant coefficient is a(d^, ep) 
with the following entries 



1 



4l> 45 

{U T Y^E) W = -j- (U 2 M E 



u D K g vlj D U* 2 ) w . 



r diag T/ r 



diag - v r * Jn M^v r \ n m) 



(50) 
(51) 



The color triplet contribution to p — > e^n accordingly reads 



T(p e+TT ) 



(V^U^M^U^u - (V UD ) u m d 



diag 



Vi} D M% a9 V?j D U* 2 )u {l + D + Ff 



where represents the VEV of the 45-dimensional representation. With the ansatz given in Eq. (47) we 
obtain the following expressions 



T{p -> e + vr°) 



m p a 



\(y UD )n\ 2 m 2 b m 2 T {l + D + Ff 



2 



|(V^jj)i2| mim 2 T (l + D + FY. 



(52) 
(53) 



These should be compared with corresponding results for the exchange of the triplet in the 5-dimensional 
representation given in Eqs. (48) and (49) to obtain 



1,2. 



(54) 



T(p -> etvr ) 45 y v v $ J 

\r oi' max x ' 

We see that the 5-dimensional triplet dominates over the 45-dimensional triplet for moderate values of U45 
if all other relevant parameters are the same. 



2. The neutral anti-lepton final state 

In order to incorporate p — > 7r + i> and p — > K + v decay modes in our study we note that one is free to sum 
over the neutrino flavors in the final state. The relevant coefficients that enter widths for these decays are 
a(d a , dp, vi) and a(d a , dp, Vi) when the exchanged state is the triplet in the 5-dimensional representation. 
To find them we need 

{D T Y z N) pi = - J- (3M^% D U;V EN + V^M^Ven)^ (55) 
where the relevant sums yield 

^(D T y%) Qi ( J D T y%)^ = 

i 

= £(3M£ ag ^[/| + VZ D U 2 Mp) aj (3M^vT D U 2 + V^M^pj (56) 
5 j 
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The upshot of these results is that widths for decays with neutral anti-lepton in the final state again depend 
only on U% as far as the mixing parameters are concerned. One can maximize amplitudes for p — > K + u 
and p — > tt + v by taking the contributions proportional to m T in Eq. (56) and using the same ansatz for U2 
as before. We find that T(p — > K + T>) dominates over widths for proton decays into charged anti-leptons. It 
reads 



2 



(57) 



This is an important result. It tells us that the most conservative bound on the scalar sector comes from the 
proton decays into neutral anti-leptons in the final state. 

If one compares the most conservative contributions of the triplets in the 5- and 45-dimensional repre- 
sentations towards p — > K + v one obtains 

Again, the 5-dimensional triplet dominates over the 45-dimensional triplet for moderate values of U45. 

We can now numerically analyze all the decay modes given in Table IX to find the current bounds on 
the triplet mass in SU (5) with symmetric Yukawa couplings. We take values of quark and lepton masses at 
Mz, as given in [44] , and neglect any running of the relevant coefficients, for simplicity. These effects can 
be accounted for in a straightforward manner. The CKM angles, when needed, are taken from Ref. [41]. 
We have randomly generated one million sets of values for nine parameters of U2 unitary matrix and five 
phases of Vud to find the bounds presented in Table X. As it turns out, it is p — > K + v that dominates in all 
instances. With that in mind we can write that the most and least conservative bounds read 

ma > 1.2 x 10- ( °) ( l0 °^-) G e V. (59) 

V0.0112GeVV V "5 I 

mA >1.5xW e—X'fm^L) Gev. ,60) 

V0.0112GeV 3 y \ v 5 J 



To summarize, if one is to maximize contributions from the triplets in the 5- and 45-dimensional rep- 
resentations towards proton decay within renormalizable SU (5) framework with symmetric mass matrices 
the current bounds on the triplet mass scale are given in Eqs. (59) and (60) if the color triplet in 5 of Higgs 
dominates in the most and least conservative scenario, respectively. In other words, any SU(5) scenario 
where the triplet scalar mass exceeds the most conservative bound of Eq. (59) is certainly safe with regard 
to the proton decay constraints on the scalar mediated proton decay. If triplet is to be lighter than that, one 
needs to explicitly check if particular implementation of Yukawa couplings allows for such scenario. The 
5-dimensional triplet dominance is determined through relations given in Eqs. (54) and (58). Finally, if the 
triplet mass is below the least conservative bound of Eq. (60) the SU (5) model is not viable. 
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Channel 


m A (GeV) 


m A (GeV) 


p -> 7r°e+ 


1.9 x 10 10 


4.9 x 10 12 


p -> 7T°/i + 


2.8 x 10 10 


5.4 x 10 12 


p -> i*T°e + 


1.7 x 10 10 


1.5 x 10 12 


p ->■ ts:V + 


2.0 x 10 10 


2.1 x 10 12 


p — > ?7e + 


1.1 x 10 10 


4.0 x 10 11 


p — > 


7.2 x 10 9 


2.7 x 10 11 


p — > 7T + J/ 


2.2 x 10 10 


8.0 x 10 12 


p — > K + v 


1.5 x 10 11 


1.2 x 10 13 



TABLE X. The least conservative (second column) and the most conservative (third column) experimental lower 
bounds on triplet mass in the 5-dimensional representation of SU (5) with symmetric Yukawa couplings. 

B. Color triplet in flipped SU(5) 

Flipped SU(5) is well-known for the so-called missing partner mechanism that naturally addresses 
scalar mediated proton decay by making the triplet scalar in the 5-dimensional representation heavy 
enough. Be that as it may, the operators associated with the triplet exchange can be suppressed with 
ease if (U T {Y W + Y WT )D) la = and {D ] c Y^*U c ) al = 0, a = 1,2, in the Majorana neutrino case. 
Note that in flipped SU (5) the Dirac mass matrix for neutrinos is proportional to the up-quark mass matrix. 
This implies that flipped SU (5) predicts Majorana nature of neutrinos. 

The operator suppression can be implemented only if the Yukawa couplings are treated as free parame- 
ters. That is not the case in the minimal realistic version of flipped SU(5) theory where Yukawa couplings 
are related to fermion masses and mixing parameters. In the minimal scenario it is sufficient to have only 
one 5-dimensional scalar representation present to generate realistic charged fermion masses. We accord- 
ingly analyze predictions of a flipped SU(5) scenario with a single color triplet state. To be able to compare 
the flipped SU (5) results with the case of ordinary SU(5) we again take Mjj^^e = ^ud E- 



1. The charged anti-lepton final state 

We get the following result for the p — > e^7r° (6 = 1,2) partial decay widths 

r(p -> e+vr ) = \(V UD Mm d - m u )\H(m 2 u + ml)\(U 2 ) 15 \ 2 (l +D + F) 2 , (61) 

where v^(= \/2 246GeV) represents the VEV of the 5-dimensional representation. Interestingly enough, if 
U2 takes the form given in Eq. (47) it would yield suppressed partial decay widths for p — > e^ir , 5 = 1,2. 
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In other words, the setup that enhances partial proton decay rates with charged anti-lepton in the final state 
in SU (5) framework suppresses corresponding rates in flipped SU (5). 

If we want to be conservative with regard to the limit on to a it is sufficient to maximize relevant decay 
widths. This can be done by taking (L^ia = 1 for p — > efir (5 = 1, 2) to obtain 

/ \ 1/2 

m A > 3.6 x 10 10 ( ~ I GeV. (62) 

V0.0112GeV 3 y 

This limit comes out to be significantly weaker with respect to the corresponding limit we presented in the 
SU(5) case. 

2. The neutral anti-lepton final state 

To find decay widths for p — > -k+v and p — > K+v channels we need to determine the form of 
a(d a , dp, Vi) and a(d a , d^ , Vi) coefficients. In the minimal model with symmetric mass matrices the rele- 
vant input reads 

( ^ (y io + y ioT )jD)ia = -^-m D M^) la , (63) 
(D T Y~ 5 N)^ = -^{V^ D M^ g U* 2 V E N) Pl , (64) 

^5 

(D^*U*)pi = --(vi D M^) pi , (65) 

^5 

with the following sum over neutrino flavors 

^(DVNUiDVN)^ = - 2 {VUM^?V^ D U. (66) 

I 

The sum over neutrino flavors in the final state eliminates dependence on U 2 — the matrix that represents 
mismatch between rotations in the quark and lepton sectors — leaving us with decay widths that depend 
only on known masses and mixing parameters. This makes minimal flipped 577(5) with symmetric mass 
matrices rather unique. We find the following limit on the triplet mass that originates from experimental 
constraints on p — > K + v channel 

/ \ l/ 2 

m A > 1.0 x 10 12 I ^ ) GeV. (67) 

V0.0112GeVV 

The fact that p — > ir + i> is also a clean channel in the sense that it depends only on the CKM angles 
and known fermion masses means that the minimal flipped SU (5) predicts ratio between T(p — > tt + u) and 
T(p K+v). We find it to be 

ffr ^ = 9.0. (68) 
T(p K+v) 
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This result is not sensitive to exact value of the nucleon matrix elements and running of relevant coeffi- 
cient. It thus represents firm prediction within the framework of the minimal flipped SU(5) with symmetric 
Yukawa couplings. 



VI. HIGHER ORDER CONTRIBUTIONS 

In the SU(5) framework the states (3,1,4/3) and A 3 E (3,3,-1/3) violate B and L and do not 
contribute to dimension-six proton decay operators at tree-level. Antisymmetry of their Yukawa couplings 
to two up quarks only allows for dimension-six operators involving c or t quarks that produce B number 
violation in charm or top decays [45], but these operators do not affect the proton stability due to large 
masses of c and t quarks. However, an additional W boson exchange opens decay channels with final states 
that are kinematically accessible to proton decay. 



A. Box mediated dimension-six operator from (3, 1, 4/3) € 45 

One possibility is to make a box diagram with a single W exchange leading to the d = 6 operator, as 
shown on Fig. 1 . In the literature, proton decay mediation involving W boson exchanges were considered 




d, s u d. s 



FIG. 1. Box diagrams with (3, 1, 4/3) state that generate d = 6 operators of flavor uudl and uddv. 

in [45^47]. We calculate the box diagram in the approximation where we neglect external momenta, how- 
ever, we keep both virtual fermions massive since the right-handed A interactions force chirality flips on 
internal fermion lines and thus the diagram would vanish if both fermions were massless. Evaluation of the 
diagrams with W and would-be Goldstones leads to gauge invariant and finite amplitude. Then we find that 
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A(3, 1,4/3) generates two effective coefficients: 



a(d a ,el 



4tt 2 



m 



ir 



Y\U ] C (Y W * -Y 10 ^)U* C 



1? L 



D ] C Y^E* C 



a(d a ,dp,Vi) 



Gf 



<i-K 2 m 2 w 4" 



m u , V ja m dk V* k J (x a , x Uj , x dk ) , 
V-^ ( yl0*_ y 10 t)[/( 



>C 



D^ C Y^E* C 



is 



pi 



(69) 



(70) 



m u , Vj a J (xa, x Uj , x £i ) 



Here, V = Vckm, while the leptonic mixing matrix has been set to unity. (For the neutrino final states one 
would need to sum over all neutrino flavors.). Mass dependence, apart from helicity flip factors, is encoded 
in function J (where x k = m^jm^) 

(y-A)y\ogy (z-A)z\ogz 



J(x,y,z) 



+ 



(y- x){y- z) (z - l)(z - y)(z - x) 

{x — 4)x log x 



(71) 



+ 



(x — l)(x — y)(x — z) 

There are two distinct regimes of dynamics in the box, depending on the presence of t quark in the loop. 
When j = 3 we expand to leading order in x £i , Xf, <C 1, and find 



J(xA,x t ,x ii ) 



1 



XA ~ X t 



XA 



xa — 1 



logXA 



x t 



1 



log x t 



When both fermions are light compared to W the J function takes the following form 



J (lA) X u . , X£.) 



1 

Xa 



XA 



log XA + 



xa - 1 x U] - x iz 



{x lt log x ti - x Uj log X Uj 



(72) 



(73) 



Contributions of the up-quark Yukawa couplings are weighted approximately by m Uj Vjd for j = 2, 3 that 
run in the box. The large mass of the t quark comes with small element Vtd that makes this product of 
the same magnitude as m c V c d- Similar cancellation between mass and CKM hierarchies occurs for the 
down-quarks and the weights obey rridVud ~ Tn s V us ~ rribV u b. 

The (3, 1,4/3) state has been identified as a suitable candidate to explain the anomalous value of the 
muon magnetic moment. One of the leptoquark couplings (D^Y^ Eg)i2 between the muon and one of 
the down quarks di must be of the order ~ 2, while other two must be < 10 -3 to suppress contribu- 
tions to other down-quark and charged lepton observables [31]. In addition, the CDF and D0 measure- 
ments of forward-backward asymmetry in ti production can be explained by a large diquark coupling 
31 ~ 2 to ut quark pair [28]. Additional couplings between uc and ct quark pairs 
are constrained by charm and top physics processes and their upper bounds are of the order 10 _1 and 10~ 2 , 
respectively [32]. Both puzzles can be explained for a mass of the leptoquark of around 400 GeV. How- 
ever, for a light mass and with the abovementioned two large couplings the proton would decay much too 
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quickly to a muon final state via dominant contribution of t quark and one of the down quarks in the box 
(c.f. Eq. (70)). Therefore, one has to find a second amplitude of equal magnitude and opposite phase in 
order to achieve cancellation between the first and the second amplitude. As explained in the preceding 
paragraph, all down quarks in the box that couple to external muon have similar weights that come from 
loop dynamics and CKM factors. As a result, the hierarchy of di contributions to the amplitude follows 
very closely hierarchy of the leptoquark couplings (D^Y^E^) a and the required cancellation cannot take 
place between the different down-quarks in the box. Likewise, cancellation between c and t quarks in the 
box diagram cannot occur for similar reason. 



B. Tree-level dimension-nine operator from (3, 1, 4/3) e 45 



The W emission from the up-type quark leads to proton decay amplitudes depicted on Fig. 2. Decays of 
this type have been already mentioned in Ref. [45]. We focus here on the final state with a single charged 
lepton whose decay width is most severely bounded experimentally. In this case the following d = 9 
effective operator is obtained 




FIG. 2. d = 6 proton decay operator induced by tree level (3, 1, 4/3) £ 45 and W exchanges. 



E 



-^GfVuoV* 



•07 



Uc (Y 10 - Y WT ^U ( 



c 



u=ct m uml 

e a bc{u^^Ld hoi )(d^aRl i )(d kl ^ iJL Lu k ) 



D ] C Y^E* C 



pi 



(74) 



Here U labels c or t quark, whereas external leptons i,j = 1,2 and down-type quarks a, f3, 7 = 1, 2, are 
all light, a, b, c, k are SU(3) color indices. R(= (1 + 7s)/2) is the right projection operator. We focus 
immediately on best constrained channels, i.e., p — > TT°£f, and we set a, /3, 7 = 1. The use of Fierz 
transformations leads to the amplitude with scalar bilinears 

8iG F 



Mi 



E 

U=c,t ' 



Uc{ Y w - Y 10T YU ( 



c 



w 



D^Y^E^ 



1; 



VudK 



ad 



(75) 



{£fRua) (d k Rd c ) (u^Ldb) 



p ) + tensor terms . 
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One can estimate the above matrix element by employing the vacuum saturation approximation. We insert 
the current dkRd c between the vacuum and tt and end up with product of pion creation and proton annihi- 
lation amplitudes. The tensor terms which are invoked by the Fierz relations in Eq. (75) cannot contribute 
in this case. The vacuum-to-pion amplitude is 



<7T° | d k Rd c | 0) 



&ck , 



(76) 



whereas the full amplitude is 



M 



p— >7T u eT 



E 



^ ™ A mu 

tabc ' 



U C (Y 10 - Y 



{ U CLd b ) {ef'Ru c 



c 



1U L 



D^Y^E^ 



, . vUdVud— — 
ii m d 



(77) 



P 



The annihilation matrix element of the proton in Eq. (77) has been most precisely evaluated using lattice 
QCD [48]. These authors have introduced operators 0^ ds = eabdu^Td^T' s c and defined constant a as 



aRu r , 



o\o%L\p) , 



(78) 



where u p is the Dirac spinor of the proton. The recent value of a obtained from lattice QCD calculation 
with domain wall fermions [43] is a = —0.0112(25) GeV 3 . The decay width is then 



r(p -> tt°( 



GlftmW X(m 2 p , ml, ml) 1 / 2 (m 2 p + m\ 



ml 



16Trm d 



rrv 



(79) 



U=c,t 



v ud [uc<y™ - y^u* c 



C\\U 



D^ C Y^E* C 



li 



rn 



where X(x, y, z) = (x + y + z) 2 — 4(xy + yz + zx). From the experimental limits in Tab. IX one obtains 
the following bounds 



[7=c,i 



CJic/ 



U=c,t 



mu 



< 2.4 x 10 



< 2.6 x 10 
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m. 



( «»<wvF cwr1 ' (80) 



-20 



(400 GeV) 2 



f— ^GeV" 1 . (81) 



A comment is in order how phenomenologically preferred values of leptoquark and diquark couplings cope 
with the above constraints. Couplings to the electrons should be small and are in particular not bounded from 



below, so the constraint from r(p — > ir°e + ) can be avoided by putting 



D ] C Y^E* C 



ii 



effectively to zero. 



On the contrary, low-energy leptoquark constraints, especially the (g — 2) M , indicate that 



D^Y^E^ 



12 
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could be large in some scenarios [31]. In this case we must require cancellation between the c and t quark 
amplitudes that occurs when 



[U C (Y W - y 10T )t?7*] 



12 



Vtd m c 
V cd m t 



2.7 x 1(T 4 x e~°- 37i . 



(82) 



[[/ c (yiO-ylOT )tC/ *] i3 

This can be achieved since the [Uc(Y 10 — Y 1QT )^U^\ 12 is only bounded from above while at the same 
time [[/c(y 10 — Y 10T )^U^\ 13 is bounded from below to satisfy observations in ti production. Finally, 
relative phase between the two couplings can be freely adjusted since it is not probed by any experimental 
observable to date. 

Finally, for the state (3, 1,-2/3) +6 present in the flipped SU(5) framework we can easily adapt the 
results obtained above since the two states are indistinguishable at low energies, provided we make the 
following substitutions 



UI(Y W * - Y^)U* C 



l r 

2 



U^Y^U, 



c 



D^ C Y^E* C 



D^ c Y l *E* c 



(83) 



To conclude, we note that despite the absence of the tree-level contribution to proton decay of the 
(3, 1, 4/3) state, weak corrections lead to proton destabilizing d = 6 and d = 9 operators. The effect of the 
d = 9 operators can be rendered adequately small even in the case of simultaneously large leptoquark and 
diquark couplings, a situation that is favored by observables in ti production and value of (g — 2)^. This is 
achieved by finely-tuned cancellation of two amplitudes. To the contrary, similar cancellation is impossible 
in the case of d = 6 operator for p — > 7r°/i + decay and we are required to suppress either all leptoquark 
couplings involving \i or all diquark couplings. We conclude that the proton decay lifetime constraint allows 
to fully address either Ap B or (g — 2)^ observable with the (3, 1, 4/3) state, but not both. 



VII. CONCLUSIONS 



Lepton and baryon number violating interactions are inherently present within grand unified theories 
and are most severely constrained by the observed proton stability. Proton decay can be mediated by vector 
or scalar leptoquarks that violate both baryon and lepton number by one unit. Vector leptoquarks that 
mediate proton decay have gauge couplings to fermions and are not readily allowed to be far below the 
unification scale. For the scalar leptoquarks, however, the freedom in Yukawa couplings gives one more 
maneuverability to realize scenarios with light scalar states. On the other hand, the very same Yukawas that 
are responsible for proton decay very often need to account for the observed fermion mass spectrum. An 
example of a setting with light leptoquark states was presented in [15, 31, 32] where the low mass of the 
state (3, 1, 4/3) had an impact on low-energy flavor phenomenology. Most notably, it was found that by 
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tuning independently the two sets of Yukawa couplings, namely the leptoquark and diquark Yukawas, one 
could reconcile the measured value of forward-backward asymmetry in tt production and the value of the 
magnetic moment of muon. 

In this work, we have classified the scalar leptoquarks present in SU (5) and flipped SU (5) grand unifi- 
cation frameworks that mediate proton decay. In both frameworks the considered leptoquark states reside in 
scalar representations of SU (5) of dimension 5, 10, 45, or 50. We integrate out the above states at tree-level 
and parameterize their contributions in terms of effective coefficients of a complete set of dimension-six 
effective operators. The mass constraint on the color triplet state contained in the 5- and 45-dimensional 
representations is then derived. The precise lower bound depends on the value of the vacuum expectation 
values of these representation. For the vacuum expectation value of 100 GeV the least (most) conservative 
lower bound on the triplet mass that originates from the p — > K + v channel is approximately 10 11 GeV 
(10 13 GeV). The corresponding bound is derived within the flipped SU (5) framework to read 10 12 GeV 
and proves to be mixing independent. Moreover, the minimal flipped SU(5) theory with symmetric mass 
matrices predicts T(p — > 7r + P)/r(p — > K + v) = 9. 

The two leptoquark states that do not contribute to proton decay at tree-level are (3,1,4/3) and 
(3,1,— 2/3) +6 in the standard and flipped SU(5) frameworks, respectively. We have estimated their 
contribution to dimension-six operators via box diagram and the tree-level contribution to dimension-nine 
operators. For the (3, 1, 4/3) state it has been found that if it is to explain both the anomalous magnetic 
moment of the muon and the tt forward-backward asymmetry, then the contribution of the dimension-six 
operator would destabilize the proton in p — > fi + n° channel. Therefore only one of the two puzzles can be 
addressed with this leptoquark state. 

Light scalar leptoquarks can be either produced in pairs or in association with SM fermions at the LHC 
and are a subject of leptoquark and diquark resonance searches [6, 49-51]. To conclude, we can expect to 
find signals of these leptoquark states at the LHC, although it seems very unlikely, in light of the constraints 
from the proton lifetime measurements, that they would be observed in a baryon number violating processes. 
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